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SUMMARY 

/ 


Aeroelastic instability phenomena of isolated open and closed rigid 
bodies of revolution free to move under elastic restraint have been 
investigated experimentally at low speeds by means of models suspended 
at zero angles of attack and yaw on slender flexible struts from a wind- 
tunnel ceiling. Three types of instability were observed - flutter simi- 
lar to classical bending-torsion flutter, divergence, and an uncoupled 
oscillatory instability which consists in nonviolent continuous or inter- 
mittent small-amplitude oscillations involving only angular deformations. 
The speeds at which this oscillatory instability starts were found to be 
as low as about one-third of the speed at flutter or divergence and to 
depend on the shape of the body, particularly that of the afterbody, and 
on the relative location of the elastic axis. 

An attempt has been made to calculate the airspeeds and, in the 
case of the oscillatory phenomena, the frequencies at which these insta- 
bilities occur by using slender -body theory for the aerodynamic forces 
on the bodies and neglecting the aerodynamic forces on the struts. How- 
ever, the agreement between the speeds and frequencies calculated in this 
manner and those actually observed has been found to be generally unsat- 
isfactory, with the exception of the frequencies of the uncoupled oscil- 
lations which could be predicted with fair accuracy. The nature of the 
observed phenomena and of the forces on bodies of revolution suggests 
that a significant improvement in the accuracy of analytical predictions 
of these aeroelastic instabilities can be had only by taking into account 
the effects of boundary- layer separation on the aerodynamic forces. 

-'-Supersedes the recently declassified NACA RM L 53 EO 7 , "An Explora- 
tory Investigation of Some Types of Aeroelastic Instability of Open and 
Closed Bodies of Revolution Mounted on Slender Struts" by S. A. Clevenson, 
E. Widmayer, Jr., and Franklin W. Diederich, 1953. 
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INTRODUCTION 


Flutter, divergence, and similar aeroelastic instability problems of 
wings and tail surfaces have been recognized for a long time. On the 
other hand, the related problem of aeroelastic instability of bodies of 
revolution (generally hereinafter referred to simply as "bodies") has 
become of interest only recently, primarily because only recently have 
external stores and fuel tanks in the shape of bodies of revolution been 
carried on high-speed airplanes, and only at high speeds do the aero- 
dynamic forces exerted on bodies at low angles of attack become suffi- 
ciently large to give rise to aeroelastic problems. 

There are several differences in the aeroelastic instability phe- 
nomena of wings and of bodies, that is, in the nature of the motions, in 
the nature of the aerodynamic forces involved, and in the nature of the 
resulting phenomena. 

The aeroelastic phenomena of wings essentially involve deformations 
of the wings themselves, whereas those of bodies are very unlikely to 
involve significant deformations of the bodies and are based, instead, 
on the deformation of the members supporting the body. For instance, a 
fuel tank carried on two struts, one behind the other, under a wing, or 
a ram jet carried similarly on supports above the fuselage can move later- 
ally as a result of the sidewise deflections of both struts in the same 
direction, and they can be yawed by a deflection of the front strut to 
one side and of the rear strut to the other. In these two degrees of 
freedom, classical flutter may occur under the proper circumstances; 
under other circumstances and involving only the yawing degree of free- 
dom, classical divergence may occur. 

The aerodynamic forces on wings at small angles of attack or under- 
going oscillations of small amplitude about zero angle of attack can 
generally be calculated with sufficient accuracy by potential-flow theory; 
they are linear functions of the angle of attack or the amplitude, respec- 
tively, and are not influenced in an essential way by the boundary layer. 
(Exceptions to this statement are the forces causing such nonclassical 
instability phenomena as stall flutter, aileron buzz, and wing buffeting.) 
The aerodynamic forces on bodies of revolution, however, are often essen- 
tially determined by the effects of viscous flow. For instance, the lift 
which is known to exist on bodies at an angle of attack in steady flow is 
due entirely to these effects, because potential-flow theory predicts zero 
lift for this case. This lift is often an intrinsically nonlinear function 
of the angle of attack. (See ref. 1, for instance.) Consequently, the 
aeroelastic instability phenomena of bodies are more likely to be of a 
nonclassical type related to stall flutter and similar phenomena than are 
the aeroelastic instability phenomena of wings. 
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Two experimental flutter investigations have been made of bodies 
mounted on wing tips (refs. 2 and 3), but no aeroelastic -instability- 
studies appear to have been made previously of essentially isolated 
bodies, that is, bodies mounted at some distance from a lifting surface 
on flexible struts which contribute no aerodynamic forces. Some instances 
where this problem arises in practice are external stores or tanks car- 
ried on struts under the wing and ram jets carried on supports on top of 
the fuselage. Also, an analysis of the aeroelastic instability of an 
isolated body may serve to shed some light on the much more complicated 
problem of aeroelastic instability of a body mounted on a wing tip. 

An investigation has therefore been conducted in order to gain some 
insight into the nature of the instability phenomena of such isolated 
bodies of revolution. A streamlined body, an open tube, and several 
bodies consisting of the tube with variously shaped end pieces were sus- 
pended from the ceiling of a wind tunnel on struts of several stiffnesses. 
The closed bodies were intended to simulate external stores or fuel tanks; 
the open tube, an unfired ram jet. In one series of tests the tube was 
also mounted on two struts covered by a fairing. The nature of the vari- 
ous types of aeroelastic instability that occurred under various condi- 
tions was observed, as were the airspeeds at which they occurred and the 
frequencies of any oscillations present. All tests were conducted at low 
speeds (Mach numbers less than 0.5) and with a range of Reynolds number 

(based on body length) from 1.5 X 10^ to 7»1 X 10^. 

In an attempt to analyze some of these results, the speed and oscil- 
latory frequency at which various types of aeroelastic instability may 
occur have been calculated by using slender -body theory for the aero- 
dynamic forces, with certain additional assumptions in the case of the 
open tube. The derivation of these forces is presented in the appendix 
of this paper, and the calculations are described therein. The calculated 
and observed results are compared and discussed. 


SYMBOLS 


a ratio of distance of elastic axis of supporting strut 

system rearward of midpoint of body to one-half length 
of body; in case of bodies consisting of tube with various 

2s l 

end pieces, the midpoint is that of tube, - 1 

L 

f frequency of oscillatory instability, cps 

f h natural lateral bending frequency of body on struts measured 

in still air, cps 
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natural frequency of yawing oscillations of body on struts 
measured in still air, cps 

structural damping coefficient (see, for example, ref. 4) 
lateral translation at elastic axis, ft 
amplitude of lateral-translation oscillation, ft 

moment of inertia of body about elastic axis of configu- 
ration, slug-ft^ 

effective lateral bending spring constant of supports with 
body mounted, lb/ft 

effective yawing spring constant, ft-lb/radian 
reduced frequency parameter, Loo/2v 

aerodynamic force per unit length along body, lb/ft 

length of body (length of tube, in case of bodies con- 
sisting of tube with end pieces), ft 

aerodynamic moment about elastic axis, ft-lb 

mass of body, slugs 

aerodynamic (lateral) force, lb 

dynamic pressure (at onset of instability, unless specified 
otherwise), lb/sq ft 

dimensionless dynamic -pressure parameter, SqVb/Ko, 
radius of body of revolution, ft 

dimensionless cross-sectional area of body, nR^L2 

coordinate along length of body, measured rearward from 
nose, ft 

distance from nose of body to elastic axis of support 
system, ft 

distance from nose of body to center of gravity of body, ft 
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V v 


volume of body, cu ft 


x a 


a 


o: 


speed (at onset of instability, unless otherwise speci- 
fied), ft/sec 

ratio of distance of center of gravity of body rearward 
of elastic axis of support system to one-half length of 


body, 


S 2 - S 1 


L/2 

yawing angle, radians 

amplitude of yawing oscillations, radians 


I dimensionless coordinate, s/L 

p density of test medium, slugs/cu ft 

P Q density of air at standard sea-level condition, 0.002378 

slugs/cu ft 

o dimensionless distance from nose of body to elastic axis, 

si/L 

to angular frequency of oscillatory instability, radians/sec 

^b. natural frequency of yawing oscillations of body on struts 

measured in still air, radians/sec 


DESCRIPTION OF TESTS 
Apparatus 


Wind tunnel .- The tests were conducted in air at variable pressures 
in the Langley 4. 5 -foot flutter research tunnel, which is of the closed- 
circular) throat single-return type. 

Models . - The two basic models, the airfoil-shaped body and the open 
tube, are shown in figure l(a). The airfoil-shaped body of revolution 
is that generated by rotating an airfoil about its chord. The ordinates 
of the airfoil are twice those of an NACA 65 -OIO airfoil and are listed 
in table 1. The open tube consisted of aluminum sheet l/l 6 inch thick 
rolled into a cylinder with a diameter of 6 inches. Various end sections 
were used in conjunction with the tube to represent closed bodies of 
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revolution with a cylindrical center section; these end sections are 
shown in figure l(b). The open tube and the airfoil-shaped body are shown 
mounted on their supports in the tunnel in figures 2(a) and 2(b), respec- 
tively. (The scales shown in fig. 2 read in inches.) 

Strut-support systems.- The models were mounted on one of four sets 
of supports which consisted of two small-diameter steel rods fixed on one 
end to a mounting plate, and on the other end to a mounting bar to which 
the models were bolted. (See figs. l(a) and 2.) The struts were designed 
to make the lateral-bending frequency of the bodies about one-half or one- 
third of the frequency of their yawing oscillation in still air. The 
natural frequency of the forward and rearward oscillations of the bodies 
was approximately six times their lateral-bending frequency. The strut 
diameters and effective spring constants of the support systems comprised 
by these struts are listed in table 2. 

The effective spring constants were obtained from the frequencies 
f ^ and fo; measured in still air and the known masses and moments of 

inertia by means of the relations 

K h = (^fh) 2 ® 


and 

% = (2rtf a ) 2 I a 


The values of K h and Kq, obtained in this way represent spring constants 
in the true sense of the word only when x^ = 0, because only then are the 

yawing and sidewise— bending oscillations uncoupled (and even then only if 
the additional apparent mass of the still air is neglected) . The values 
given in table 2 are averages of the values obtained with different bodies 
for Xq, = 0 (except for the values listed for struts A, for which fre- 
quency measurements were made for Xq = 0.10 and 0.14-) . 

For one series of tests with the open tube, the struts B were covered 
with a fairing of aluminum alloy 1/32 inch thick, which extended about 
1 inch ahead of the front strut and 1^ inch behind the rear strut and was 

attached to the struts along their entire lengths but was not attached to 
the mounting plate or mounting bar. The airfoil obtained in this manner 
was about 3/16 inch thick at and between the two struts, had a rounded 
leading edge, and a sharp trailing edge. 
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Strain-gage instrumentation .- The only instrumentation, apart from 
the usual instrumentation required to measure the tunnel speed and den- 
sity, consisted of electrical-resistance strain gages mounted on the 
roots of the struts in such a way as to measure strains due to sidewise 
deflections of the struts. The output signals of the gages on the front 
and the rear struts were amplified separately and fed to two channels of 
a multiple recording oscillograph and also to an oscilloscope for immedi- 
ate visual observation. 


Tests 

General procedure .- The procedure for each test was to increase the 
tunnel speed slowly and at the same time the angle of yaw of the model was 
adjusted (by yawing the mounting plate by means of a mechanism outside 
the tunnel) to give zero lift and moment on the body, as indicated by the 
strain-gage outputs. When some type of instability occurred, the strain- 
gage outputs were recorded, the type of instability was noted, the tun- 
nel conditions were observed, and the test was discontinued, except in 
some instances when the instability was not violent and it was desired to 
study it further. 

The model-strut combinations tested in this manner are listed in 
the left half of table 3 ; also listed are the model mass, model moment 
of inertia, elastic-axis and center-of-gravity location, as well as the 
measured still-air frequencies of each configuration. The tests are 
divided into several series for the sake of convenience in referring to 
them. 


Tests on the streamlined body .- In series I, the streamlined body 
was mounted on struts A, the most flexible ones, and the tests were con- 
ducted at various air densities. In series II, the same body was mounted 
on the somewhat stiffer struts B. Only the nonviolent uncoupled oscil- 
latory instability occurred even at the highest air- density used in these 
tests (substantially sea-level density). The tunnel speed was increased 
successively to several values beyond that at which this instability first 
occurred, the air density being kept substantially constant at the sea- 
level value. Test series II consists of measurements of the frequencies 
of the oscillations at these air speeds. . The streamlined body was also 
mounted on the still stiffer struts C, but no aeroelastic instability of 
any kind occurred; consequently, this experiment is not listed in table 3 * 

Tests on the open tube with miscellaneous end sections .- The open 
tube with various end sections was mounted on struts B, C, and D, that 
is, on all but the most flexible struts. Series III consists of the tests 
made with the various configurations at constant air density. Series IV 
consists of tests made at various densities by using the tube with hemi- 
spherical sections at both ends mounted on struts C. 
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Tests on the open tube on unfaired struts . - Tests made at various 
air densities on the open tube mounted on the most flexible struts com- 
prise series V, and tests made at constant air densities on the open tube 
with various center-of-gravity positions mounted on the next stiffer 
struts comprise series VT. The center of gravity was varied by attaching 
narrow bands of lead l/l6 inch thick to the inner surface of the tube, so 
that the mass, the moment of inertia, and, hence, the still-air frequencies 
were changed as well. The open tube was also mounted on struts D, but no 
instability occurred at any speed up to and including the speed at which 
this test was discontinued, namely, 536 feet per second. This test is, 
therefore, not listed in table 3* 

Tests on the open tube on faired struts .- The tests which constitute 
series VII are those made at constant air density on the open tube mounted 
on struts B with the fairing attached. The location of the center of 
gravity was varied in the same manner as in series VI. 


RESULTS 

Presentation of Results 


The results of the tests are presented in the right half of table 3 
and some of these results are presented in figures 3 to 9* 

The speeds listed in table 3 for tests which led to flutter or diver- 
gence are those at which these instabilities first occurred. Similarly 
the frequencies listed for the tests which led to flutter are those at ’ 
the flutter condition. The air speeds listed for the tests of series IV 
are those observed when the instability first occurred, and the frequen- 
cies are those observed at that speed. Similarly, in the tests of 
series II, the first speed listed is that at which the instability first 
occurred, but the other speeds are merely speeds above the first at which 
the frequency of the oscillation was measured. The last speed is that at 
which the tests of series II were discontinued; the nature of the insta- 
bility phenomenon did not appear to change in the speed range covered. 

The flutter-speed coefficient 2 v/Iao ( -j, the dimensionless dynamic 
pressure at flutter q*, and the frequency ratio f/f a pertaining to the 
tests of series I are plotted in figure 3 against the density ratio p/p Q . 

The frequencies measured in series II are shown plotted as a function of 
airspeed in figure L. In figure 5 are shown the speed and dynamic pres- 
sure (both in dimensionless form) at the onset of the yawing oscillations 
observed in the tests of series IV as functions of the density ratio p/p Q 

The flutter-speed coefficients 2 v/La^, and the frequency ratios f/f a 

pertaining to the tests of series V are shown plotted as functions of the 
density ratio P/P Q in figure 6. 


NACA TN 3308 


9 


In figure 7 the speed coefficients for which flutter or divergence 
occurred are shown plotted against Xq, for the tests of series VI. The 

points representing the various tests are not connected because the body 
mass and moment of inertia (and, hence, the still-air frequencies) were 
not constant in these tests. Figure 8 consists of a similar plot made 
for the tests of series VII. 


Observed Flutter and Divergence Phenomena 

Several types of aeroelastic instability were encountered. In the 
tests of series I and V, as well as in some of those of series VI and VII, 
coupled flutter similar to classical bending-torsion flutter was encoun- 
tered, except that "bending” and "torsion" were lateral motion and yawing, 
respectively, in the case of these bodies. These two types of motion 
were distinguished by observing the strain-gage outputs. If the outputs 
of the front and rear gages had been in phase and of the same magnitude, 
the motion would have been from side to side only, without yawing, but 
actually this type of motion was not observed in the tests. When the 
gage outputs were l 80 ° out of phase and of the same magnitude, the motion 
was a pure yawing oscillation, and, when they were out of phase by any 
other angle, the motion was a coupled lateral-motion and yawing oscillation 

In one test of series VTI a combined flutter and divergence insta- 
bility was observed, not unlike the type of phenomenon which a wing may 
experience if its flutter occurs in a mode which involves very little 
bending; the tube tended to diverge to the stops after a few oscillations 
of increasing amplitude. As in all tests where divergence was observed, 
when the body began to diverge in yaw it moved over laterally as well as 
under the action of the side forces brought into play by the yawed attitude 

The flutter frequencies were ill-defined occasionally, particularly 
when the body at its flutter condition was also close- to a divergence con- 
dition; that is, flutter then occurred so suddenly that no definite fre- 
quency could be obtained from the strain-gage record. 


Observed Uncoupled Oscillations 

An instability was observed in the tests of series II, III, and IV. 
This phenomenon consisted in continuous or intermittent, self-excited, 
small- amplitude yawing oscillations, usually with fairly well defined 
frequency. When the oscillations were intermittent they started up at 
random intervals rather than subsiding and increasing in a regular fashion, 
as do oscillations with beats. In two of the tests in which such oscil- 
lations occurred the frequency was insufficiently defined to be measured. 
This phenomenon differed from flutter not only in the fact that it involved 
small constant amplitudes but also in the fact that, unlike flutter, it 
involved no bending deformations of the struts and, hence, no lateral 
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motions of the bodies. For lack of a better name this phenomenon is listed 
as "yawing oscillations" in table 3 and will be referred to as such here- 
inafter . If the body were mounted in such a way that the struts were hori- 
zontal, as would be the case if a body were mounted on the side of a fuse- 
lage, this phenomenon would consist of angle-of-attack or pitching 
oscillations . 

Some of the tests of series III which resulted in yawing oscillations 
were continued to speeds above that at which the oscillations started; the 
first speed listed is then the one at which the instability first occurred, 
and the second is that at which the tests were discontinued. No change 
was noted in the nature of the phenomenon within this range of speeds. 

The two frequencies listed for these tests correspond to these two speeds; 
the values of p, q, and q* are those which correspond to the first 
speed. In the tests of series III in which no instability occurred, and 
also in the first test of series VI, the speed listed is that at which the 
tests were discontinued. 


Results of Calculations 

Some of the results of the instability calculations described in the 
appendix are also listed in table 3 and are shown in figures 4, 6, J, 
and 8 For series I and II the calculated flutter speed was infinite 
that is, the calculations did not predict flutter for any finite speed. 
Nor was it possible to calculate a finite speed at which self- sustained 
yawing oscillations could exist, but the frequencies of the yawing oscil- 
lations of the body in a moving airstream could be calculated and are 
shown in table 3 and figure 4. 


Flutter speeds and frequencies were calculated for the tests of 
series V; the results are shown in table 3 and figure 6. Flutter speeds 
and frequencies as well as divergence speeds were also calculated for the 
tests of series VI and VII and are listed in table 3; the speed coeffi- 
cients are also shown in figures J and 8. 


DISCUSSION 

A Note Concerning the Speed and Dynamic -Pressure Parameters 


Two dimensionless parameters have been used in order to compare the 
results of the various tests. The first of these is substantially the one 
commonly used in flutter work, the flutter parameter or speed coeffi- 
cient v/bo where b is the half -chord and is here replaced by one- 
half the (basic) body length, so that the parameter becomes 2v/ln) a . It 
involves a measure of the dynamic pressure (in the term v, although the 
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air density is not taken into account), as well as of the structural and 
inertia properties involved in angular deformations (combined in the 
still-air yawing frequency a) a ) . The other parameter used herein is one 

often used in static aeroelastic work and may be referred to as the 
divergence parameter or dimensionless dynamic pressure; the form of this 

2qV h 

parameter used in this paper is q* = ■ ■ . It represents the ratio of 

% 

the aerodynamic moment per unit angular displacement to the elastic 
restoring moment per unit displacement and involves implicitly the dynamic 
pressure as well as the structural and aerodynamic properties pertinent 
to angular deformation, because the factor 2 is the value of the moment 
coefficient (referred to the volume of the body) per unit angular dis- 
placement according to thin-body theory. This theory applies only to 
closed bodies, but, if it is extended to open bodies on the basis of the 
assumption that the flow inside the body has the same velocity as that 
outside and that the rear end of the body acts like a sharp trailing edge, 
the same value is obtained for the moment coefficient because the effects 
of these two assumptions cancel each other. 

The main advantage of the flutter parameter is that it includes some 
dynamic or inertia effects; on the other hand, the advantage of the 
divergence parameter is that by virtue of its explicit inclusion of aero- 
dynamic properties it serves as a more precise definition of certain 
instability phenomena. For instance, flutter and divergence can occur 
over a wide range of values of the flutter parameter, but, although flut- 
ter may occur at almost any value of the divergence parameter, divergence 
should occur at values of this parameter near unity, (if 2 is the correct 
value of the moment coefficient per unit angular displacement or if the 
correct value is used in the definition of the divergence parameter 
instead of 2, divergence will occur when the parameter is 1.) Also, as 
shown in the appendix, the divergence parameter appears to play an impor- 
tant part in determining the frequency of the uncoupled yawing oscillations. 

Both parameters, therefore, have some advantages and, in view of the 
exploratory nature of the investigation, both have been used in attempts 
to analyze the observed phenomena. 


Flutter and Divergence 

Flutter and divergence similar to the classical instability phenomena 
on wings were both observed on the aerodynamically clean bodies, the 
streamlined body and the open tube with and without fairing on the struts, 
as may have been expected, because under certain conditions these bodies 
have linear aerodynamic forces which may be expected to give rise to 
phenomena similar to the classical instability phenomena of the wings. 

The streamlined body .- The streamlined body fluttered in the tests 
of series I, in which the density was varied but all other parameters 
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held constant, at speeds which corresponded to a wide range of the flutter 
parameter 2v/In^, hut to values of q* which varied only between O.79 

and 0.94. The flutter frequency was substantially constant in these tests. 
(See fig. 3.) 

The quantity q* varied relatively little in these tests because 
it contains the air density, which was the only variable in these tests. 

The small variation of q* in the tests can be considered to be due to 
the change in the mass ratio, that is, the ratio of the body mass to the 
mass of the displaced air. Inasmuch as the body would have diverged at 
q* = 1 if the actual moment on this body in steady flow were that pre- 
dicted by thin-body theory, flutter occurred in these tests at dynamic 
pressures from 6 to 21 percent lower than the theoretical dynamic pres- 
sure at divergence. 

When attempts were made to calculate the flutter speed of the stream- 
lined body, the aerodynamic forces predicted by thin-body theory were 
found to be incapable of predicting a finite flutter speed. For bodies 
with a fineness ratio of about the moment due to angle of attack and 
normal force due to steady rotation predicted by thin-body theory are 
about 25 and 4-0 percent higher, respectively, than those predicted by 
exact potential-flow theory. The values predicted by potential-flow 
theory are, in turn, somewhat higher than the actual values as the result 
of boundary- layer separation. Inasmuch as quantitative errors in the 
predicted forces would tend, by themselves, to result only in an incorrect 
flutter speed, the fact that the predicted flutter speed does not even 
exist suggests that aerodynamic forces must be involved which are not 
predicted by this theory. Such forces are the normal force due to angle 
of attack and the moment due to steady rotation. These forces are zero 
according to thin-body and exact potential-flow theory, but actually they 
do exist; often they vary linearly with angle of attack and rate of rota- 
tion, respectively, and thus represent the type of forces required to 
cause classical flutter. 

In the light of this discussion, prediction of the flutter speed of 
bodies of revolution thus requires a knowledge of the effects of the 
boundary layer and of the phenomena associated with its separation on 
the aerodynamic forces. Hence, the main shortcoming of thin-body theory, 
insofar as the prediction of flutter is concerned, consists in the inability 
of potential-flow theory in general to predict some of the critical forces 
involved in these phenomena rather than in the degree to which thin-body 
theory approximates exact potential-flow theory. In divergence, however, 
only the moment due to angle of attack is involved; therefore, the correct 
magnitude of this force is important. Inasmuch as the moment predicted 
by thin-body theory is about 30 or 40 percent higher than the actual 
value, the dynamic pressure at divergence would be that much higher than 
that estimated on the basis of thin-body theory; therefore, in the tests 
of series I, the highest speed reached probably corresponds to about 
80 percent of the true divergence speed. 
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The open tube on unfaired struts .- In the tests of series V the open 
tube fluttered at all densities. The values of 2v/Lo^ decreased from 
10.32 to 5*30 with increasing density, and the values of q* appeared 
to increase gradually from O.516 to 0.579* except for the second test of 
the series for which q* was somewhat higher than expected from the 
other three tests. (See table 3*) These trends are the same as those 
observed in the case of the streamlined body. In contrast to the stream- 
lined body, however, flutter speeds could be calculated for the open tube 
on the basis of thin-body theory (with the additional assumptions of 
unobstructed flow through the tube and of the flow continuing in the same 
direction after leaving the tail end of the tube rather than reassuming 
the free-stream direction). These speeds are in fairly good agreement 
with the measured ones. 

The observed frequencies did not vary with the density, a fact which 
was also noted in the case of the streamlined body. The calculated fre- 
quencies did not vary with density either but were about 30 percent higher 
than the observed ones. 

On the basis of these comparisons it appears that the aerodynamic 
forces are taken into account in a qualitatively correct manner, but that 
quantitatively they must be improved considerably before they can be used 
to predict flutter speeds and frequencies correctly for ducted bodies. 

In the tests of series VI, the tube fluttered or diverged in all 
cases, except in the first test, in which the speed of the test was not 
carried to a high enough value. Flutter occurred when the center of 
gravity was at or behind the elastic axis, and divergence, when it was 
at or in front of the elastic axis (see fig. 7)* when it was at the 
elastic axis, the tube fluttered at the higher mass and diverged at the 
lower mass . This trend agrees with the trends noted in the tests of 
series I and V, because in these tests there appeared to be a tendency 
to approach divergence as the air density increased and, hence, the mass 
parameter decreased. The values of 2v/ln^ at instability varied from 

5. 47 to 7-95 and those of q* between 0.^59 and 0.75^- (see table 3)* 

The values of q* at divergence tended to be higher than those at flut- 
ter, as may be expected in view of the nature of q* as, primarily, a 
divergence criterion. No such distinction can be made in the case of 
the values of 2v/Lo^, at instability, both the highest and lowest values 

of which correspond to divergence. 

The calculated speeds at instability are in fairly good agreement 
with the measured speeds for the rearward locations of the center of 
gravity, but the flutter speeds predicted for forward center -of -gravity 
locations are much too low. The two measured flutter frequencies are 
substantially below the values calculated for those two cases. The cal- 
culated divergence speeds are consistently higher than the measured ones 
by about 20 percent on the average; therefore the moment coefficient per 
unit angular displacement (the only aerodynamic parameter which enters 
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the calculation of the divergence speed) must he about 40 percent higher 
than that estimated by means of the modified thin-body theory used herein. 
Actually the tube under consideration is not sufficiently slender to 
justify the use of thin-body theory, and its abrupt changes in cross- 
section at the nose and tail violate certain assumptions inherent in 
thin-body theory; also, the validity of the two additional assumptions 
is doubtful. Hence, the inability of the modified thin-body theory to 
predict quantitatively useful results is not surprising. 

A more accurate potential-flow solution for the flow through and 
about an open tube would be difficult to obtain, and its validity would 
still be open to question because the flow inside the tube would be 
significantly affected by the boundary layer in the inside walls. The 
available solution for a ring airfoil is inapplicable to this case 
because the tube is far too slender for this theory. The most promising 
solution, therefore, appears to be the use of a semiempirical method for 
estimating the aerodynamic forces required in flutter analysis. Such a 
method might consist in retaining thin-body theory but modifying the two 
additional assumptions, that is, by estimating the magnitude of the 
forces which are, in effect, concentrated at the rear of the afterbody 
and the extent to which the flow decelerates inside the tubes on the 
basis of measurements of the moment and normal force due to angular dis- 
placement. In divergence calculations only the moment per unit angular 
displacement is required, of course. 

The open tube on faired struts . - The aeroelastic instability phenom- 
ena of the tube were substantially unchanged by the addition of the 
fairing; apparently, the increase in the aerodynamic forces was canceled 
in effect by the increase in the stiffness of the configuration. Flutter 
still tended to occur at the further rearward position of the center of 
gravity, and divergence, at the forward positions. (See table 3 a nd 

Q # ) The agreement between the calculated and measured speeds was 
poor, the calculated values being much too low, and the calculated fre- 
quency corresponding to the one measured frequency was also much too low. 
The first two tests of this series serve to illustrate the difficulty of 
estimating the speed at which aeroelastic -instability phenomena occur; 
under identical test conditions, the model diverged in one case at 
322 feet per second and exhibited some symptoms of flutter at 481 feet 
per second in the other case, flutter apparently having been suppressed 
the first time. 

The fact that the calculated values of the divergence speed are 
lower than the observed values indicates that, inasmuch as the aerodynamic 
moment on the tube due to angular displacement was probably underesti- 
mated, as in the tests of series VI, the forces on the fairing were 
greatly overestimated, as might be expected to be the case because unmodi- 
fied two-dimensional theory was used to estimate these forces. The exact 
calculation of the mutual interference effects of the tube and the wing 
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tip represented by the fairing represents a very difficult problem for 
which no solution has been found to date. For the time being, therefore, 
approximations similar to those made herein must be used, although they 
could probably be improved by resorting to empiricisms based on some 
measured results . 


Yawing Oscillations 

The streamlined body; speed at which oscillations start .- When in 
the tests of series II the streamlined body was tested on stiffer struts 
and with an elastic-axis location 5 percent of the body length further 
forward than in the tests of series 1 , no flutter was observed} instead, 
the body experienced the self -excited yawing oscillations described pre- 
viously. These oscillations started at a speed which is relatively low 
compared to those at which flutter occurred in the tests of series I; 

it corresponds to = 3.76 and q* = 0 . 211 , and the speed at which 

Ho 

a 2 v 

the tests were discontinued corresponds to = 6. JO and q* = 0.630, 

whereas at the same density flutter would have occurred in the tests of 
series I at values of 2.vjlxa^ and q* of about 8.0 and 0.9* respec- 
tively. Therefore, if the tests of series II had been continued to a 
speed some 15 or 20 percent higher than that at which they were discon- 
tinued, flutter would probably have occurred if the values of 2v/liu Q , 
and q* at flutter are assumed not to differ much between the two test 
series. The fact that oscillations occurred in the tests of series II 
but not in those of series I is probably the result of the difference in 
elastic-axis locations, as will be shown later . 

In order to calculate the speed at which yawing oscillations might 
start, an attempt was made to solve the equations of motion for the case 
of zero lateral displacement. For this case there are two differential 
equations with one unknown function. With the aerodynamic forces given 
by thin-body theory, however, the equations can have a solution only when 
the airspeed is zero. Therefore, the nature of the air forces must again 
be different from that assumed to yield equations which are compatible 
at nonzero airspeed. 

From a physical point of view, self -excited oscillations can occur 
only if the net damping is zero. In the equation for the equilibrium 
of the moments as set up in the appendix, there is no damping term, 
because the aerodynamic damping moment according to thin-body theory is 
zero and the structural damping was assumed to be zero. If, however, 
the aerodynamic moment were actually small and negative, at a certain 
speed it would be just large enough to counteract the structural damping, 
and self-excited oscillations would start at this speed. At higher 
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speeds, the net damping would be negative, but nonlinearities in the 
aerodynamic forces might prevent the oscillations from diverging. Inas- 
much as these aerodynamic forces would then be due entirely to deviation 
from potential flow, they would probably vary with Reynolds number so 
that the question whether and when a body might experience self -excited 
oscillations probably depends on the Reynolds number involved, as well 
as on the structural damping. 

The streamlined body; frequency of oscillations .- Although the equa- 
tions set up for the analysis of the yawing oscillations do not furnish 
a solution for the speed at the onset of the oscillations they do give 
an indication of the frequency of the resulting oscillations, that is, 
of the "natural" frequency in moving air. This frequency is expressible 
as a product of the still-air yawing frequency, which involves the dynamic 
characteristics of the system, and a correction factor which, except for 
a generally negligible dependence on the mass ratio, is a function only 
of the static aeroelastic characteristics represented in the param- 
eter q*, namely vr - q*. That the trend of the calculated frequencies 
agrees well with the trend of the frequencies measured for the stream- 
lined body may be seen from figure 4. However, the rate of decrease 
with speed is less than predicted. This fact suggests that, inasmuch 
as the expression for the frequency is obtained directly from the equa- 
tion representing the equilibrium of the moments on the body, the esti- 
mated aerodynamic moments on the streamlined body are too large, which 
is true, as previously noted. 

Miscellaneous bluff bodies; effect of elastic -axis location on the 
speed at which oscillations start .- Yawing oscillations occurred in the 
majority of the tests of series III. The speeds at which they started 
correspond to values of 2v / /Lj^ ) , and q* much lower than those at which 

the streamlined body fluttered in the tests of series I and of about the 
same magnitude as those at which the oscillations of the streamlined body 
started in series II. (See table 3.) Some of the tests of series III, 
for instance, the last-listed one, resulted in no instability at what 
appear to be fairly large values of 2v/Iu^ and q*; however, all values 

of 2v/Lo> a listed in table 3 are based on the length of the basic tube. 

If they had been based on the actual lengths of the various bodies, 

2v/lcx a would have been about 4 for the last-listed test, which value is 
about the same as that at which the oscillations of the streamlined body 
began (3.7 6) and much lower than that at which that body could have flut- 
tered at this density (about 8). Also, although the value of q* (O.76O) 
at which this test was discontinued is high compared with that at which 
the oscillations of the streamlined body started (0.211), it is lower 
than that at which that body would probably have fluttered at this density 
(about 0.9). The values of q* in table 3 for series III are based on 
the actual volume, as are those for the other series. 
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The speeds at which oscillations started, in the form of the parame- 
ter 2v/La^, are shown plotted in figure 9( a ) a s functions of the dimen- 
sionless elastic -axis location a. The importance of the elastic -axis 
location in determining the speed at which the oscillations start may he 
deduced from the equation set up ‘in the appendix for the equilibrium of 
moments in the yawing-oscillations phenomenon, which contains the dimen- 
sionless elastic-axis location a. 

In figure 9(a) the upward -pointing arrows on several sketches refer 
to cases in which no instability was observed; therefore, any instability 
would have had to occur at values of 2v/Lo) a greater than those shown, 

which represent those at which these tests were discontinued. Figure 9(a) 
also contains one point representing the conditions at which the stream- 
lined body started to oscillate in the tests of series II. The sketches 
in figure 9(a) all imply an airstream approaching from the left. 

The only body for which sufficient information was available to 
deduce the effect of elastic -axis location on the speed at which the 
oscillations started is the one with hemispherical nose and tail. In 
figure 9( a ) this body is represented by four points which appear to lie 
on a smooth curve which has a minimum at a = 0. For the body with hemi- 
spherical nose and streamlined tail, two points were available. These 
points are connected by a curve based on the pattern exhibited by the 
body with hemispherical nose and tail which represents the probable 
variation of the speed coefficient with elastic-axis location, although 
the minimum of the curve may not be at a = 0 as shown. For each of 
three other bodies (including the body used in the tests of series I 
and II), one point was available which represented the onset of oscilla- 
tions and one point which indicated only that the oscillations, if any, 
would have to start at higher speeds except that, for the steamlined 
body, the second point represented the flutter condition. Curves repre- 
senting the estimated variation of speed at the onset of oscillations 
with elastic -axis location are shown for two of these bodies as well. 

These approximate curves indicate that the speed at which oscillations 
start is lowest when the elastic axis is near the midpoint of the body. 

By using the estimated variations as a guide, that is, by estimating 
what the speed would have been if the elastic -axis -location parameter a 
had been zero, the effect of the body shape on the speed at which oscil- 
lations tend to start can be divorced from that of the elastic -axis loca- 
tion. 


Miscellaneous bluff bodies; effect of body shape on the speed at 
which oscillations start.- Inasmuch as a part of the large differences 
in the speed coefficients shown in figure 9(a) for the various bodies is 
due to the fact that the coefficients were based on the length of the 
basic tube, the speed coefficient will be considered to be based on the 
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actual length of the body in the following discussion. The various 
bodies then fall into several main classes. Oscillations appear to start 
at the lowest speeds when the tail of the body is hemispherical, regard- 
less of whether the nose is hemispherical, square (cover plate), or a 
small cone (for all of which 2v/Lo) Q , is about 1.3 )> or a streamlined 

body (for which 2v /Lo^, is about 1.5). The class with the next highest 

range of speeds at which oscillations start is that with tails consisting 

of the streamlined shape or the large cone. If for this group the nose 

is hemispherical or square, 2v^In> Q , is about 2.5; if the entire body is 

streamlined, 2V/LCDQ, is about 3; and if the nose of the body is a pointed 
streamlined shape or a large cone, 2v / /Lo^ is greater than 4. The class 

for which oscillations start at the highest speeds is, surprisingly, that 
with a square tail (cover plate); if for this group the nose is hemispheri- 
cal or square, 2v/Lx> a is about 3 to 4, and, if the nose is a pointed 

streamlined shape, 2v^Ia> Q , is greater than 4. All these numerical values 

must be used with caution, of course. 

In general, then, the speed at which oscillations start appears to 
depend to a large extent on the shape of the tail of the body, a hemi- 
spherical tail being the least favorable in that it oscillates at the 

lowest speeds, a streamlined tail or large cone used as a tail being much 
better, and a square tail being most favorable in this respect. The 
shape of the nose has almost no effect when the tail is hemispherical but 
has some effect in the other cases, a square or hemispherical nose being 
worst, a conventionally streamlined (rounded) nose being better, and a 
pointed streamlined nose, best. The aerodynamically cleanest configura- 
tion, the streamlined body, occupies a relatively favorable place; the 
speed at which it may start oscillating can be increased further by 
replacing the rounded streamlined nose by a pointed one. 

The only way in which the numerical values given for 2v/lxa a in the 

preceding discussion can be related to those corresponding to the classi- 
cal instability phenomena is by noting that the streamlined body fluttered 
at values of 2v/Li> q; of about 8. As a rule of thumb, then, based on these 

very limited data, an aerodynamically clean body may be expected to start 
oscillating at speeds as low as about one-third its flutter or divergence 
speed. To attempt a similar correlation for bluff bodies would be futile, 
because these bodies do not experience the linear aerodynamic forces on 
which the classical instability phenomena are premised. 

The reasons for the relative behavior of the various bodies are as 
yet somewhat obscure. The effects of the nose on the speed at which 
oscillations start are probably associated with separation at the nose, 
because the less disturbance caused by the nose, the more favorable the 
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conf iguration. Similarly, the behavior of the tail can probably be 
explained by the effects of separation at or ahead of the tail. The 
unfavorable effects of the hemispherical tail, for instance, are probably 
the results of the rather large projected area of the tail (projected on 
a plane through the axis of the body), which is exposed to the separated 
flow. The relatively favorable square tail, on the other hand, has no 
projected area exposed to this flow. The streamlined tail probably 
causes relatively little separation and thus is relatively favorable if 
the nose is favorably shaped; this effect corresponds to the observed 
effects . 


In view of the complicated nature of these phenomena there appears 
to be little hope of arriving at theoretical methods of predicting the 
speed at which the oscillations start, at least in the case of bluff 
bodies which are unlikely to be exposed to high-speed airstreams anyway; 
in the case of more or less streamlined bodies, empirical approaches may 
prove fruitful. In either case the oscillations do not appear to be vio- 
lent, and, if they are undesirable, they can always be eliminated by 
stiffening the supports and in many cases merely by changing the body 
shape or the elastic -axis location, and possibly also by using vanes 
mounted on the body. 


Miscellaneous bluff bodies; frequency of oscillations .- The ratio 
of the frequency of the oscillations to those in still air is shown 
plotted in figure 9(b) as a function of the dimensionless dynamic pres- 
sure q* corresponding to the speed at which the oscillations started. 
The arrows in figure 9(b) refer to cases in which no frequency was meas- 
ured and serve merely to call attention to the fact that oscillations 
did start on the given model at the indicated value of q*. As in fig- 
ure 9( a ), one point in figure 9(b) represents the conditions at which 
the streamlined body started to oscillate in the tests of series II. 
Also, as in figure 9( a )> the sketches in figure 9(b) imply an airstream 
approaching from the left. 


The frequencies shown in figure 9(b) agree fairly well with the theo- 
retical curve (which neglects the effect of the mass ratio). In view of 
the fact that many of the assumptions made in the analysis are violated 
by the bluff bodies, this agreement is better than may have been expected. 
All points in figure 9(b) pertain to tests at densities close to sea 
level. The results of the tests of series IV at various densities follow 
a similar pattern, although the range of q* covered by these tests is 
small. (See fig. 5 .) The frequency of the oscillations can thus be 


estimated on the basis of the relation 


f - jrr 

■‘■a 


q* with fairly good 


accuracy. If the body is aerodynamically clean and the speed relatively 
high, the accuracy of this formula can be improved by replacing the fac- 
tor of 2 in the definition of q* by a better value of the moment coef- 
ficient per unit angular displacement. At low speeds, the frequency may 
be expected to be substantially the same as the still-air frequency. 
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CONCLUSIONS 


1. Streamlined bodies and open tubes mounted on thin flexible struts 
which do not contribute any aerodynamic forces have been found to diverge 
and flutter; flutter tended to occur for relatively far rearward locations 
of the center of gravity and for relatively high mass ratios (body mass to 
mass of displaced fluid). 

2. Flutter could not be predicted for the streamlined body by using 
aerodynamic forces based on potential theory. For the open tube, the 
calculated flutter and divergence speeds did not agree well with the 
measured values; the discrepancies are believed to be due to the intrinsic 
shortcomings of potential-flow theory. The analysis of unsteady aero- 
elastic effects of bodies of revolution therefore appears to require a 
knowledge of the boundary-layer and separation effects on the unsteady 
forces on these bodies. 

3. Closed strut-mounted bodies of revolution appear to be subject 
to a nonclassical instability which consists in self-excited nonviolent 
oscillations which appear to start, in the case of aerodynamically clean 
bodies at least, at speeds about one -third that at flutter or divergence 
for the given body. 

k. The speed at which the oscillations start for a given body depends 
on the elastic -axis location, being lowest when the elastic axis is 
located near the midpoint of the body. This speed is also determined to 
a large extent by the shape of the body, particularly of the tail end. 

For bodies with hemispherical tails, the oscillations start at low speeds 
but, for bodies with streamlined and, particularly, with squarely cut off 
tails, they start at relatively high speeds; the optimum nose shape appears 
to be a pointed streamlined shape and the next best, a conventional stream- 
lined nose shape. 

5 . The mechanism which causes these oscillations is as yet unknown, 
although negative aerodynamic damping moments appear to be required. The 
speed at which the oscillations start cannot, therefore, be predicted at 
present; its calculation apparently requires a knowledge of boundary- 
layer and separation effects on the unsteady aerodynamic forces. The 
frequency, however, can be estimated from a simple formula involving the 
frequency of the oscillations in still air and the ratio of the given 
dynamic pressure to that at divergence. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 12, 1953. 
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APPENDIX 

CALCULATION OF AEROELASTIC INSTABILITIES 


Equations of Motion 


The equations of motion of a rigid body mounted on flexible supports 
and performing combined yawing and lateral oscillations are 


-IojO: - m ^2 - s-j^'h - (l + ig)KQO! + = 0 


( 1 ) 


-m^S2 - s-j^ - mh - (l + ig)K h h + P = 


( 2 ) 


where the dots designate differentiation with respect to time, and where 
P is the aerodynamic force and Mq, the aerodynamic moment (positive in 

the direction of positive h and a, respectively) . The manner in which 
P and Mq, were calculated for the bodies considered in this paper is 

described in the following sections. 


Aerodynamic Forces on Closed Bodies 

The aerodynamic forces on closed bodies of revolution performing 
unsteady motions in supersonic flow have been calculated by linearized 
potential-flow theory (see refs. 5 and 6). The aerodynamic forces on 
bodies of revolution in steady incompressible flow can easily be calcu- 
lated by potential-flow theory (by using sources and sinks, for instance, 
as in ref. 7); according to this theory, the normal force is zero. The 
exact calculation of the aerodynamic force and moment for unsteady motion 
by potential-flow theory, however, is quite difficult both at subsonic 
and supersonic speeds, and in view of the fact that they are known to be 
influenced to a large extent by the effects of viscosity, a large expendi- 
ture of effort in calculating them is hardly warranted. In the absence 
of any means of taking the effects of viscosity into account for unsteady 
motions, a simple approximation to potential-flow theory, namely slender - 
body theory, has therefore been used for the purpose of the calculations 
described herein. (See ref. 8, for instance, for an outline of a slender - 
body theory in quasi-steady flow.) 

The assumption made in slender -body theory is that the momentum of 
the flow in a plane perpendicular to the free stream is the same as it 
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would be if this flow were two-dimensional. This assumption implies 
that the derivative of the radius with respect to distance along the 
length of the body is small (which implies, in turn, that the body is 
very slender, that is, the body has a high fineness ratio), and, also, 
that the angle of attack and any motions are small. (For a fuller dis- 
cussion of these assumptions from the mathematical point of view, see 
ref. 9 .) 

The momentum of the flow about a circular cylinder for a unit length 
along the cylinder is equal to the product of the rate of motion of the 
cylinder and the apparent mass, which is equal to the mass of air dis- 
placed by a unit length of the cylinder, or p«R2. At any section of a 
body of revolution the rate of motion relative to the component of the 
free-stream velocity normal to the axis of the body is va + h + (s - sp)( 

this value is within the approximations implied in slender -body theory. 
Therefore, the momentum per unit length along the body of the flow per- 
pendicular to the body at this section is pjtR2(s) jva + h + ^s - s-^j Ctj . 

The force exerted by the body on the fluid per unit length along the 
body is equal to the time rate of change of the momentum per unit length, 
the rate of change being that along the path of a particle, that is, 

— . But within the approximation implied in slender -body theory. 


Dt 



Dt 8 t 8 x 


_ 8 _ 

8 x 


Therefore, the force per unit length along the body is 


l 



( 3 ) 


Hence, if the body is performing oscillations defined by 


icut 


( 4 ) 


icot 



( 5 ) 
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then 


c 


7 /• x 8 \ 
l = -p [ lO) + V — 1 


= -pj[2i 


a 1 


«R 2 (s) 


va + iaii 0 + ico^s - S]_) 


O' 


icot 


^2iawrtR‘ 


'(s) - go 2 (s - s 1 )«R 2 (s) + v 2 — itR 2 (s) + 
' 8s 


rvWs - s, nR 2 (s) 

^ 1 ks 


a o + 


|-ob 2 R 2 (s)rt + ivcjo — jtR 2 (s)lh 

L 8s J 0 


icut ( C \ 
e (6) 


The aerodynamic force P (positive in the same direction as h) 
and moment ^ (positive in the same direction as a) are then 


P = / l ds 

Jo 


= -pVu < 


iow + (jo 


( S 1 " S ) 


a_ - arh 
0 o 


iajt 


(7) 


and 


%c=J (s - 8l )l ds 
u 0 


= P V v 


1 L ' 


v 2 + cjo 2 (s 2 - 2 S 2 S + Sj 2 ) 


^o + 


nr 


(s - sj + i 


ivgo 


V e 


ioot 


( 8 ) 


where 


s = — / sjtR 2 (s) ds 

Vb^O 
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and 


s 2 = — f 1 s 2 «R 2 (s) ds 


b o 0 


or, in dimensionless form, 


p V 2 L 2e i art 


2 iklQ + 4 k 2 (alQ - Ip) 


a - 2k 2 I Q ^2. 
o u L 


(9) 


and 


M 


a 


2 t 5 iajt 
pv L^e 


I 0 + 4k 2 (a 2 ^ - 2al 1 + I 2 ) 


a o + 


iklQ + 2k 2 (ip - al 0 ) 


2h, 


( 10 ) 


where I Q = f S(|) d|; I = f |S(|) d|; and I p = f | 2 S(|) d| . 

J 0 Jo Jo 


For the airfoil-shaped body of revolution, these three values are 
I 0 = 0.01626, Ip = 0.00674, and I 2 = 0.00335- 

For steady flow (co = 0) equations (7) and (8) give the known results 
of slender -body theory: 


P = 0 

(11) 

% = 

(12) 


Aerodynamic Forces on Open Bodies 

In attempting to calculate the forces on the open tube in a similar 
manner, several problems arise. If the recovery factor is assumed to be 
100 percent, that is, if the velocity of the flow in the tube is assumed 
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to be equal in magnitude to the free-stream velocity, the combined appar- 
ent mass of both the internal and the external flow at a given section 
of the body is just twice the apparent mass of the external flow used 
heretofore. However, one of the assumptions of slender -body theory, 
namely, the one concerning continuity of the radius along the length of 
the body, is violated at the nose and tail section of the tube. Conse- 
quently, as a result of the abrupt changes in cross-sectional area, not 
only the concentrated forces predicted by slender -body theory at the nose 
and tail section but also the distributed forces predicted on the remain- 
der of the tube are open to question. 

For lack of a better theory a modified slender-body theory has been 
used to calculate aerodynamic forces on the open tube. The modification 
consisted in disregarding the concentrated forces on the tail section on 
the premise that both the external and the internal flows leave the 
trailing edge of the tube tangentially and are not realined with the free 
stream. This assumption is equivalent to the Kutta condition of subsonic 
airfoil theory and is used also in the application of slender-body approxi- 
mation to airfoil theory. (See ref. 10.) The assumption is thus, essen- 
tially, that the exit section of the tube acts like the trailing edge of a 
wing of very low aspect ratio. 

With these approximations equations (6), ( 9 ), (10), (ll), a ^d (12) 
become, in the case of an open cylindrical tube. 



where 6(s) is the delta function which represents the concentrated 
loads, and 


r 





( 14 ) 
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< 2ik(l - a 2 ) - 0 - Lk 2 Q- - o + a 2 ) 


a + 
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-2k 2 f i - a I - Iko 
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(15) 


P oj= 0 = -2prtR 2 v 2 Q! 


( 16 ) 


M. 


. = 2pitR 2 s-,v 2 a 

u a >=0 


(IT) 


Calculation of Flutter Speeds and Frequencies 

For the streamlined body, the force and moment coefficients given 
by equations ( 9 ) and (10) were introduced into equations (l) and (2) 
with a and h given by equations (k) and ( 5 ), and an attempt was made 
to solve these equations by the conventional methods of two-degree-of - 
freedom flutter analysis. (See ref. 4, for instance.) However, no 
solution was found to exist; therefore, if the aerodynamic forces were 
correct, the body would not experience flutter at any finite speed. For 
the open tube the force and moment coefficients given by equations (l4) 
and ( 15 ) were substituted into equations (l) and (2). In this case 
flutter speeds and frequencies did exist and the computed values are 
given in table 3 and are shown in figures 6 to 8. 

For the tube on the faired struts the forces and moments given by 
equations (l^) and ( 15 ) were used for the tube proper. For the fairing 
the force and moment were assumed to be given by two-dimensional theory 
at any section and were obtained from reference 4. The aerodynamic inter- 
action between tube and fairing was thus neglected, as was the effect of 
the finite span on the forces experienced by the fairing. A Rayleigh-Ritz 
type of analysis was used with two modes, a linear torsion and a parabolic 
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bending mode. These modes were selected on the basis of the consider- 
at ion that the struts deflect somewhat as shown in the two following 
schematic front views: 



Lateral displacement 



The flutter speeds and frequencies calculated in this manner are given 
in table 3 and are also shown in figure 8. 


Calculation of the Frequency of the Uncoupled Yawing Oscillations 


In order to determine what characteristics, if any, of the yawing- 
oscillation type of instability could be predicted, cognizance was taken 
of the fact that these oscillations did not involve any bending deflec- 
tions; hence, h was set equal to zero in equations (l), (2), ( 9 ), 
and (10) and, for the sake of convenience, the structural damping coef- 
ficient was assumed to be zero as well. The following equations resulted 


r 


2 2 3 

V* + pv L 

Iq + 4k 2 (<?Iq - 2oI 1 + I 2 ^ ► 




( 18 ) 
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and 
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A solution of equation (l8) can be written as 


(19) 



1 


pl5 



( 20 ) 


This solution furnishes no information concerning a speed at which the 
oscillations may start but, instead, gives the frequency at which the 
body will tend to oscillate if it is yawed away from its equilibrium 
position and then released. This frequency is also the dominant fre- 
quency of the response of the body to random excitation. The term 

pL fo 2 1 - 2crl + in the denominator of the expression is the ratio 

I \ 0 1 2/ 

a 

of the moment of inertia of air at free-stream density contained within 
the body to the moment of inertia of the body alone, both taken about the 
elastic axis. This term is inversely proportional to the mass ratio; it 
depends on the shape of the body to some extent but is substantially 
independent of the elastic-axis location. Except at extremely low mass 
ratios, this term is negligibly small; for the airfoil-shaped body for 

instance, it is 0.0020 — and 0.0019 — for the elastic-axis location 

p o P 0 

used in the tests of series I and II, respectively. 

The numerator of the expression under the radical in equation (20) 
is equal to 1 - q*, where q* is defined by 


2qV b 


( 21 ) 
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Therefore, if the small term in the denominator is neglected, equa- 
tion (20) can be written as 


a 


= )Ji - q* 


(22) 


In order to determine the speed at which these oscillations should 
occur equation (l8) must be solved simultaneously with equation (19) . 
However, the only real solution of equation (19) is v = 0 and, in 

addition, ( aI () ~ ^l)* ^ erefore > if the aerodynamic forces 

given by equations (9) and (10) were correct, oscillations could occur 
only at zero airspeed. These oscillations would then be the ordinary 
still-air yawing oscillations, the condition on x a being that necessary 

to uncouple the yawing from the sidewise-bending mode. 


Calculation of Divergence Speed 

Inasmuch as divergence is a static instability phenomenon, the speed 
required to diverge the body can be found by setting h, go, and q 
equal to zero in equations (l) and (10) or, more simply, setting go = 0 
in equation (l8) . Thus, 


or 


^divergence 


I K 

a_ 

phh 


0 



divergence 


K 

= _QL 
2 V, 


(23) 


(24) 


The parameter q* defined in equation (21) is thus equal to the ratio 
of q to divergence* ^he dynamic pressure at divergence being that 

calculated by using slender -body theory. 


The divergence speeds for the tests of series VI and VII are given 
in table 3 and are shown in figures 7 and 8. For series I and II the 
value of divergence ma y ’ be obtained from equation (24), and for any 
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other test it can be obtained from the values of q* and q given in 
table 3* Therefore, 


^divergence - for series I and II 

= 29A for series I 
= 1-92 for series II 

= — , in general 
q* 
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TABLE 1.- ORDINATES USED TO GENERATE 
THE AIRFOIL-SHAPED BODY 


s/L 

r/l 

O.OO5 

0.0154 

• 0075 

.0186 

.0125 

.0237 

.025 

.0315 

.05 

.0435 

• 075 

.0529 

.10 

.0608 

.15 

.0732 

.20 

.0829 

•25 

.0901 

• 30 

.0954 

•35 

.0985 

.4o 

.0999 


.0993 

.50 

.0960 

•55 

.0906 

.60 

.0829 

.65 

.0736 

.70 

.0631 

• 75 

.0516 

.80 

.0397 

.85 

.0277 

• 90 

.0162 

• 95 

.0061 

1.00 

0 


R 
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TABLE 2.- DIAMETERS AND EFFECTIVE SPRING 
CONSTANTS OF STRUTS 


Strut 

Strut diameter, 
in. 

K h 

% 

A 

O.OLO 

10 

20 

B 

.100 

50 

25 O 

C 

.166 

110 

430 

D 

.251 

200 

820 


TABLE 5.- RESULTS OF AEROELASTIC -INSTABILITY TESTS 


4 =- 


Body 

Nose Tail 

Strut 

m 


a 


f h 


Type of 
instability 

p x 1(>3 

V 

exp 

v anal 

^xp 

f exp 

^anal 


<kl 

q exp 









Test serit 

is I 










0.1301 


0.0690 


-0.28 


2.82 


0.54 

.83 

1.13 

1.50 
1.90 

2.50 


34i 

277 

235 

207 

187 

172 


31 A 
31.8 
31.2 
32.2 
33-2 
37-0 


1.70 
1-75 

1.71 
1.69 
1.60 
(a) 


15 A 

12.5 

10.6 
9-3 
8.4 
7.8 


0.798 
.808 
• 793 
.818 
.834 
.940 


Test series II 


B 

0.1300 

0.0670 

-0.18 

0 

4.00 

10.1 


Yawing oscillation 


2.33 

2.31 

2.28 

2.26 

2.21 

2.19 


299 

340 

387 

443 

505 

532 


104 

133 

171 

222 

282 

310 


9.25 

8.95 

8.77 

8.47 

8.06 

7-75 


8.99 

8.65 

8.18 

7.51 

6.64 

6.18 


3.76 

4.29 
4.88 
5.38 

6.30 
6.70 


0.210 

.269 

.345 

.448 

.569 

.626 


Test series III 


B 

0.1378 

0.1014 

-0.44 

0.44 

3.08 

7-9 

C 

.1378 

.1034 

.24 

-.24 

4.03 

10.2 

D 

•1535 

.1091 

-.20 

.20 

5-85 

14.4 

C 

.1463 

.1264 

0 

.06 

4.18 

9.2 

D 

.1463 

.1274 

-.20 

.26 

5.61 

12.9 

C 

.1635 

.1626 

0 

.24 

3.96 

8.0 

C 

.1635 

.1910 

.44 

-.20 

4.11 

7-5 

c 

.1392 

.0904 

0 

-.04 

4.34 

10.9 

D 

.1392 

.1045 

-.20 

.16 

5.85 

14.4 

D 

.1383 

.0926 

-.20 

.20 

5-91 

14.4 

D 

.1394 

.0950 

-.20 

.16 

5.86 

14.7 

D 

.1348 

.0926 

-.20 

.20 

5.87 

14.9 

C 

.1677 

.1630 

.28 

0 

3-92 

8.1 

D 

.1383 

.0940 

-.20 

.20 

5-90 

14.8 

D 

.1436 

.1180 

-.20 

• 30 

5.65 

13.2 

C 

.1677 

.1738 

0 

-.26 

4.00 

8.0 

C 

.1635 

.1766 

0 

-.24 

3.92 

7.8 

C 

.1776 

.2100 

0 

-.18 

3.16 

7.2 


Yawing oscillation 
Yawing oscillation 
Yawing oscillation 

Yawing oscillation 
None 

Yawing oscillation^ 
Yawing oscillation b 

Yawing oscillation b 
None 


None 

Yawing oscillation 
None 

Yawing oscillation b 
Yawing oscillation 
None 
None 

Yawing oscillation 


2.28 

2.36 

2.31 

2.32 
2.22 

2.22 

2.21 

2.05 

2.06 


2.30 

2.24 

2.04 

2.31 

2.27 

2.21 

2.18 

2.28 

2.15 


234 

156-233 

223-269 

235-301 

479 

245-324 

435 

312-412 

556 

390 

223-256 

524 

211-301 

223-256 

368 

424 

156-256 

479 


62 

29 

57 

64 

254 

67 

209 

109 

296 

175 

56 

280 

51 

56 

150 

196 

28 

247 


7-5 

10 . 3 - 10.1 

14.1-13-9 

8.58-8.34 


6.15 

(a) 


(a) 


14.7 

7-8 

14.7-14.5 
7. 5-7* 4 


3.84 

1.95-2.91 

1.97-2.38 

3.25-4.16 

4.73 

3.90-5.16 

7.39 

3.65-4.82 

4.75 

3.45 

1.93-2.22 

4.48 

3.32-4.73 

1.92-2.20 

3-55 

6.75 

2.55-4.18 

8.48 


0.276 

.075 

.077 

.181 

.377 

.190 

.593 

.266 

.378 

.229 
.072 
• 335 
.137 
.073 
.215 
.526 
•079 
.760 


a Not obtained. 
b Intermittent . 
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TABLE 3.- RESULTS OF AEROELASTIC -INSTABILITY TESTS - Concluded 


Body 

Nose Tail 

Strut 

m 

la 

a 

*a 

f h 

f a 

Type of 
instability 

p x 105 

v exp 

v anal 

^exp 

f 

exp 

f anal 

<&) 

' “ 'exp 

Wanal 

* 

^ exp 








Test series 

IV 


















- 

0.23 

423 



20.6 

10.5 



5.07 


0.053 











.42 

314 

— 

20.7 

10-7 


3.78 


.054 











.61 

233 

— 

1 6.6 

10.5 


2.80 


.043 











.80 

194 

— 

15.1 

10.7 

— 

2.34 


• 039 



0.1578 




4.35 

10.6 



.99 

1.19 

1.38 

171 

160 

160 

— 

14.5 

15.2 

17.7 

10.4 

10.7 

10.7 

— 

2.05 

1.92 

1.92 


.037 
• 039 
.046 

(] ' ' D 

C 

0.1301 

0 

0 

Yawing oscillations < 

















1.57 

162 

— 

20.6 

10.7 

— 

1.94 


.053 











1.77 

152 

— 

20.5 

10.6 


1.82 


.053 











1.96 

158 

— 

24.4 

10.6 

— 

1.90 


.064 











2.15 

145 

— 

22.6 

10.6 

— 

1.74 


.058 











2.31 

129 

— 

19.2 

10.4 

— 

1-55 


.050 








Test series 

V 


















' 

0.55 

195 

140 

10.5 

1.51 

1.99 

10.32 

7.42 

0.516 



0.1410 

0.0864 


0.14 


2.40 



1.12 

147 

103 

12.1 

1.54 

2.00 

7.80 

5.50 

• 594 


A 

0 

1-37 

Flutter ^ 


1.55 

117 

88 

10.6 

1.58 

1.96 

6.19 

4.68 

.521 











2.35 

100 

71 

11.8 

1.56 

1.99 

5.30 

3.76 

• 579 








Test series 

VI 











0.1237 

0.0681 

0 

0 

3.92 

10.1 

None 

2.31 

323 



121 




2.04 


0.471 



.1410 

.0875 

0 

.08 

3.10 

8.4 

Flutter 

2.14 

377 

c 296 

152 

(a) 

6.8 

5.71 

4.48 

.596 



.1410 

.0708 

0 

0 

3.H 

9.2 

Divergence 

2.14 

396 

C 312 

168 


d 7- 3 

5.47 

6.75 

.659 



.1410 

.0957 

0 

.14 

3.07 

8.2 

Flutter 

2.14 

369 

c 304 

146 

5.50 

6.7 

5-73 

4.72 

•575 


B 

.1580 

.1217 

0 

.26 

3.01 

7.4 

Flutter 

2.14 

353 

c 323 

133 

(a) 

6.1 

6.06 

5.56 

.522 



.1750 

.1475 

0 

.36 

2.81 

7.0 

Flutter 

2.14 

331 

c 35l 

117 

4.60 

5.6 

6.00 

6.38 

.459 



.1580 

.1217 

0 

0 

3.00 

7.3 

Flutter 

2.14 

389 

c 276 

162 

(a) 

6.2 

6.78 

4.82 

.636 



.1750 

.1475 

0 

-.14 

2.83 

6.8 

Divergence 

2.14 

424 

c 177 

192 


*6.6 

7.95 

9.13 

.754 



.1580 

.1217 

0 

-.26 

2.96 

7.5 

Divergence 

2.14 

375 

c 124 

150 

— 

*6.0 

6.37 

8.28 

.589 


Test series VII 




0.1410 

0.0708 

0 

0 

5.41 

13.7 

Divergence 

2.14 

522 

284 

292 



4.86 

2.64 

1.146 



.1410 

.0708 

0 

0 

5.41 

13.7 

Flutter and divergence 

2.14 

481 

£135 

248 


<14.7 

4.48 

61. 25 

• 974 



.1410 

.0875 

0 

.08 

5.30 

13.0 

Flutter 

2.14 

486 

*113 

253 

(a) 

3-9 

4.51 

1.11 

• 993 


U 

.1580 

.1217 

0 

-.26 

3.40 

11.1 

Divergence 

2.14 

427 

f 221 

195 


^5.3 

4.90 

3.26 

.766 



•1750 

.1475 

0 

-.14 

2.97 

8.1 

Divergence 

2.14 

393 

?99 

165 


<*2.1 

6.16 

4.46 

.648 



• 1750 

.1475 

0 

.36 

2.87 

8.1 

Flutter 

2.14 

332 

*6l 

118 

4.30 

2.0 

5.20 

• 96 

.463 


a Not obtained. ^"vNACA 

^Calculated flutter speed; the calculated divergence speed is 487. 

^Calculated flutter frequency. 

^Struts covered with fairing simulating thin airfoil. 

^Calculated flutter speed; the calculated divergence speed is 284. 
g Based on calculated flutter speed. 
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Hemisphere section 


ON 



Elastic axis 



Cover plate section 



(a) Airfoil- shaped body, basic 
open tube, and supports. 


(b) Bodies of revolution used as nose 
and tail sections for the basic 
tube . 


Figure 1.- Dimensions of the various bodies and the support system. 
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(a) Open tube. 



(b) Airfoil- shaped body. 

Figure 2.- Models mounted in the Langley 4. 5-foot flutter research tunnel. 
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Figure 3-- Plot of experimental speed, dynamic-pressure, and frequency 
parameters at flutter for the airfoil-shaped body. (Table 3, test 
series I.) 
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Figure 4.- Plot of frequency of yawing oscillations against airspeed for 
the airfoil- shaped body. (Table 3, test series II.) 
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Density ratio, — 

Po 

Figure 5-- Plot of speed coefficient and dimensionless dynamic pressure 
against the density ratio for the body with hemispherical nose and 
tail. (Table 3, test series IV.) 
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Figure 6.- Plot of experimental and analytical speed and frequency 
parameters against the density ratio for the open tube. (Table 3, 
test series V.) 
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Figure 7*- Plot of speed coefficients at flutter and divergence against 
dimensionless center-of-gravity position for the open tube on unfaired 
charts. (Table 3, test series VI.) 
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Figure 8.- Plot of speed coefficients at flutter and divergence against 
dimensionless center-of-gravity position for the open tube on faired 
charts. (Table 3, test series VII.) 
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(a) Speed coefficients at onset of (t>) Frequency ratios at onset of 

yawing oscillations. yawing oscillations. 


Figure 9-- Speed coefficient and frequency ratio at onset of yawing 
oscillations of various bodies of revolution. (Table 3, test 
series III.) 
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